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ABSTRACT 

We report on radiative hydrodynamic simulations of moderate and strong solar flares. The flares 
were simulated by calculating the atmospheric response to a beam of non-thermal electrons injected 
at the apex of a one-dimensional closed coronal loop, and include heating from thermal soft X- 
- - - ray, extreme ultraviolet and ultraviolet (XEUV) emission. The equations of radiative transfer and 

I statistical equilibrium were treated in non-LTE and solved for numerous transitions of hydrogen, 

helium, and Ca ii allowing the calculation of detailed line profiles and continuum emission. This 
work improves upon previous simulations by incorporating more realistic non-thermal electron beam 
models and includes a more rigorous model of thermal XEUV heating. We find XEUV backwarming 
contributes less than 10% of the heating, even in strong flares. The simulations show elevated coronal 
and transition region densities resulting in dramatic increases in line and continuum emission in both 
the UV and optical regions. The optical continuum reaches a peak increase of several percent which is 
consistent with enhancements observed in solar white light flares. For a moderate flare (~M-class), the 
dynamics are characterized by a long gentle phase of near balance between flare heating and radiative 
cooling, followed by an explosive phase with beam heating dominating over cooling and characterized 
^ . by strong hydrodynamic waves. For a strong flare (^X-class), the gentle phase is much shorter, and 

' we speculate that for even stronger flares the gentle phase may be essentially non-existent. During the 

\ explosive phase, synthetic proflles for lines formed in the upper chromosphere and transition region 

CO ' show blue shifts corresponding to a plasma velocity of ^^120 km s~^, and lines formed in the lower 

' chromosphere show red shifts of ~40 km s~^. 

, Subject headings: methods: numerical — radiative transfer — Sun: atmosphere — Sun: flares 

o 

^ ■ 1. INTRODUCTION 

P ,■ 

I ■ During a solar flare, large quantities of energy are transferred between the corona and chromosphere through thermal 
O conduction, non-thermal particle beams, radiation transport and mass motions. Any attempt to accurately model the 
^ optical and UV radiation produced during a flare must incorporate these methods of energy transfer. There has been 
considerable effort put into modeling the atmospheric response to heating from non-thermal particle beams. .Brownl 
' ((l,973) and Emslic ( 197^ derived an analytic expression describing electr on energy deposition through Coulomb colli- 
sions . This expression has be e n employed in a numb er of flare simulati ons llRicchiazzi fc Canfieldlll983HMacNeice et alJ 
T98l iNagai & Emsliej [1981 iFisher. Canfleld. ^A /IcClvmont 198^ iPeres et all Il987t iMariska. Emslie. fc Lil 119891: 
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Abbett fc Hawlevl ll999D . Understanding how energy is transferred by mass motions and ra- 



diative transfer requires solving the equations of radiation hydrodynamics. Many hydrodynamical simulations have 
focused on the upper atmosphere whe re emission is optically th i n, and radiative transfer is rel a tively straightforward . 
Notable in vestigations of this type are lCheng eiTal] l|198ail984|) : IMariska. Emslie. fc Lil l| 19891) : IPeres fc R.ealel l|1993D : 
lR,eale fc"P ercs (1995). However, we seek to model optical emission which originates in the chromosphere where radia- 
tion is neither optically thin nor in LTE with the surrounding plasma. Therefore, we must employ the full non-LTE 
formulation of radiati ve transfer. 

The simulations of iFisher. Canfield. fc McClvmonil l|198fi|) were among the flrst to model the radiative hydrody- 
namic response in a coronal loop to a beam of non-thermal electrons. Howev er, their models were li mited in that 
radiative transfer was solved using an escape probability method. Th e models of | Hawlev fc Fisher! l)1994j) incorporated 
non-thermal beam heating into the radiative transfer code, MULTI (!Carlsson"198d). This had the advantage of effi- 
ciently solving the equation of radiative transfer, but did not include detailed hydrodynamics, and t herefore could not 
simula te the effects of mass motions on the energy balance or on the line profiles. The models of I Abbett fc Hawlevi 
l)1999t) (hereafter, AH99) incorporated both hydrodynamics and radiative transfer, and are therefore well suited for 
understanding the optical emission produced during flares. In this paper we extend the models of AH99 by more 
accurately modeling the XEUV radiation and electron beam heating. 

In section 2 we describe the method by which the radiative hydrodynamic equations are solved. We also describe the 
procedure by which we generate the initial atmosphere and our new methods for incorporating soft X-ray and electron 
beam heating into the simulations. In section 3 we discuss the dynamics of our simulations and present detailed line 
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and continuum profiles. The limitations of these models and future improvements are described in section 4, and in 
section 5 we present our conclusions. 

2. METHOD OF SOLUTION 

The basic method of solution is described in detail in AH99 and lAbbefr3 ljl998|) . Here we briefly outline the procedure. 
The equations of hydrodyn amics, population con servation and radiative transf er are solved simu l taneously on aone- 
dimensional adaptive grid l|Dorfi k Drurvlll987j) using the RADYN code of iCarlsson fc SteinI 1)1994 119951 Il997|) . 
The code has been modified to include the effects of flare heating from a non-thermal electron beam, heating from 
high temperature soft X-ray, extreme ultraviolet and ultraviolet emission (hereafter, we refer to these collectively as 
XEUV emission) and optically thin cooling resulting from brcmsstrahlung and coUisionally excited metal transitions. 
Hydrodynamic effects due to gravity, thermal conduction and compressional viscosity are included as described in 
AH99. We have additionally modified the code to include the double power-la w electron beam energy distributions 
recently observed in solar flares with the RHESSI satellite ijHolman et al.ll2b03D . We in corporate the effects of XEUV 
heating from a large number of high temperat ure lines using results from the CHIANTI ijDere et al.lll997HYoung et alJ 
I2003|) and ATOMDB databases ( S mith et alH oOl'). 

We obtain the non-LTE solution to the population equations for a six level with continuum hydrogen atom, nine level 
with continuum helium atom and six level with continuum Ca ii atom using the technique of ISch armcr fc CarlssonI 
()1985(1 . The bound-bound and bound- free transitions which are calculated in detail are listed in Tables [5] All lines 
are calculated assuming complete redistribution, but for the Lyman transitions the effects of partial redistribution 
are approximated by truncating the line proflles at 10 Doppler widths. Other transitions are included as background 
opacity in LTE. The opacities are computed from the Uppsala package of Gustaf sson (1973 ). We include optically thin 
radiative cooling from bremsstrahlung and coronal metal transitions using emissivities from the ATOMDB database. 
The equations are solved on an adaptive grid using 191 grid points in depth, 5 in angle and up to 100 in frequency 
for each transition. We have found it necessary to employ an adaptive grid in order to resolve large shocks which 
propagate through the atmosphere. 

2.1. The Preflare Atmosphere 

A magnetic flux tube is typically modeled as a semi-circular loop with footpoints embedded in the photosphere 
and apex in the corona. By symmetry we need only consider one leg of the loop which we take to be a cylinder of 
constant cross section. In this approximation the atmospheric height, z, is the only spatial degree of freedom. We 
use the PF2 model of AH99 as the preflare atmosphe re. The PF2 atm osphere was generated by adding a transition 
region and corona to the model atmosphere of Carlsso ^ fc SteinI l|1997|) . To resolve the steep gradients present in the 
transition region PF2 required 191 grid points. Constant non-radiative quiescent heating was applied to grid zones with 
photospheric column mass (i.e. column mass greater than 7.6 g cm^^) to balance the energy losses in the photosphere. 
The upper boundary was set at 10"* km and held at 10^ K. With these boundary conditions and no external sources 
of heating, the atmosphere was allowed to relax to a state of hydrodynamic equilibrium. 

2.2. Electron Beam Heating 

In the standard flare reconnection model, electrons located near the loop apex are accelerated to high energies and 
travel downward through the atmosphere depositing energy and producing hard X-ray bremsstrahlung radiation. For 
the simulations presented here, we assume that the electron beam is the primary source of flare heating in the lower 
solar atmosphere. Recently, Holman ct al. ( 2003) analyzed hard X-ray spectra observed with RHESSI to obtain the 
injected electron energy spectrum at 20 s intervals throughout the 23 July 2002 X-class solar flare. The injected 
electrons are found to have a double power law energy distribution of the form 

p 



to{Eo) - -. — r- < ; \<S„-5, / s-5„ 

En{5u~2)-(^) i5^-5A{[^) (f) iorEo>Es 



where T is the electron energy flux that enters the magnetic loop; E^ is the cutoff energy below which the X-ray emission 
is assumed to be thermal; and Eb is the break energy where the distribution shifts from spectral index 5i to 5,,,. Using 
Equa tion ID we calculated the energy deposition rate as a function of column depth using the method of lEmsha l)1978l 
[T98l|) which assumes that the energetic particles are slowed and deposit their heat through Coulomb interactions. T o 
account for the changing partial ionization of hydrogen we have adopted the treatment of iHawlev h Fisheil l)1994(l . 
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The resulting heating rate as a function of column depth is, 
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•ylN) = xA + (1 — x)A', X is the hydrogen ionization fraction, A and A' are the Coulomb logarithms 
defined in jilmslic ( 1978J, and /io is the cosine of the pitch angle of the beam; /io = 1 in our simulations. The quantity 
Bx^{S/2, 1/3) is the incomplete beta function for Xc — N/Nc where Nc — iiqE^ /3Kj is the maximum column depth 
penetration of a beam with cutoff energy Ec- is defined similarly for Eb- N*{N) is the equivalent column depth 

in a fully ionized plasma defined by N*{N) = j{N')/ AdN' and N* is the cutoff depth in a fully ionized plasma. 

The parameters Si, Su, F , Ec and Eb are obt ained by fitting the h ard X-ray spectra observed during a flare to models 
for thick target bremsstrahlung as described in lHolman et alJ 1)20031) . In a future paper, we will use the time-dependent 
beam parameters they have obtained to simulate the 23 July 2002 solar flare. In this paper, we take 6i, Ec, and 
Eb to be constant with values of 3.0, 4.0, 37 keV and 105 keV, respectively, corresponding to the peak of the 23 July 
2002 flare. 

Figure ^ shows the initial energy deposition rate of the electron beam and compares it with the initial deposition 
rate from AII99. AII99 used a single power law electron energy spectrum with a spectral index, 6 = 5, and a cutoff 
energy, Ec = 20 keV. The lower spectral indices and the higher energy cutoff we employ result in energy deposition 
that extends deeper in the atmosphere and has a broader spatial extent. 

Non-thermal electrons affect the atmosphere not only by energy deposition but also throug h direct collisional ion- 
ization . We estimate collisional ionization due to non-thermal electrons using the technique of iRicchiazzi &: Canfieldl 
IpSl . They find the non-thermal collisional ionization rate from the ground state to be 

C„t = 3.78 X lO^T-^^^T (3) 

(1 - x)nH 

where uh is the hydrogen number density and Qi"'' is the beam energy deposition rate due to collisions with neutral 

hydrogen and is related to the total energy deposition rate defined in Equation|21by qI""* = Qe(l — x)A' /j. They also 
find that ionizations from the ground state dominate, so we neglect non-thermal collisional ionizations from excited 
states. 

2.3. Soft X-Ray, EUV and UV (XEUV) Heating 

Mass motions during fiares raise the density in the transition region and corona which can greatly increase the 
emissivity of XEUV photons from these regions. The outward directed photons we detect as elevated emission, while the 
downward directed photons cause "backwarming" : heating of the lower atmosphere due to increased photoionizations. 
We have used the ATOMDB database to determine the thermal volume monochromatic emissivity as a function of 
wavelength and temperature. We include emissivities calculated for approximately 34,000 transitions at 37 temperature 
points ranging from 10** K to 10^ K. The transitions have been gathered into 14 wavelength bins ranging from 1 A 
to 2500 A. Table El lists the range and emissivity-weighted central wavelengths of each bin. When calculating these 
emissivities, we have taken care not to include transitions that are already solved in detail in RADYN (cf. Table 

Since the XEUV emission produced in low-temperature, optically-thick layers is negligible, we assume that all such 
emission originates from the higher temperature, optically-thin layers. The XEUV flux incident on a layer with optical 
depth Ti, is given by, 

S^{U)E{t^ - T^) dU -2tt J SME{t^ - U) dU , (4) 

where Si, is the source function, given by the ratio of the emissivity to the linear extinction coefficient. For layers close 
to the layer under consideration (r < d/VS, where r is the separation between layers and d is the loop cross-sectional 
diameter), E in the above integral is E2, the second exponential integral. However, in order to avoid overestimating 
the emission from distant layers (r > d/VS) we follow ICan fc Fangi l|1990D : iHawlev fc Fisheil l|1994D : and AII99 in 
treating the distant layers as point sources. In that case E is given by, 

where /i' is defined in Equation 10 of AII99. 
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The XEUV photons heat the lower atmosphere through ph otoionization. The av e rage s cale height of our model 
atmosphere is greater than the thermalization depth scale of iHenoiix fc Nakagawa I ijlQTTT . so we assume that the 
energy of the photoionized electrons is transformed entirely to heat in each atmospheric layer. The volumetric heating 
rate is given by, 



where rij is the number density for ion j, Xj is the ionization potential, aji, is the photoionization cross section, and 
the sum is over approximately 150 ions. The cross sections and number densities were obtained from the CHIANTI 
database. 

Our method for calculating the XEUV heating is similar to the soft X-ray heating of AH99, but differs in the extent 
of the atomic transitions included. AH99 considered emission only between 1 A and 250 A divided into seven bins, 
while we include 14 bins in the range 1 A to 2500 A, i.e. including substantial EUV and UV heating. For example, 
C III A977 A, which was not included in AH99, is a significant contributor in our backwarming calculations. Figure|21 
compares our XEUV heating rate with that of AH99 at a comparable time late in a moderate flare. Our XEUV heating 
rate is approximately 10 times larger than the rate in AH99. 



We carried out two flare simulations, corresponding to moderate beam heating with an electron flux, T = 10^° ergs 
cm^^ s^^ (the FIO flare) and strong beam heating, with T = 10^^ ergs cm~^ s~^ (the Fll flare). The FIO flare was 
evolved for 226 s and the Fll flare for 15 s. Note that these simulations are comparable to the FIO and Fll flares 
in AH99 which evolved for 70 s and 4 s, respectively. In the following sections we discuss the dynamics and emission 
resulting from these two flare simulations. 



Figure 13 shows the evolution of the temperature and density stratification along with the electron beam energy 
deposition rate, the ionization fraction and the electron density during the FIO flare. Initially the beam penetrates 
into the chromosphere to a depth of 0.89 Mm above the solar surface (defined as z = at T5000 — !)■ The upper 
chromosphere rapidly heats but plateaus at a temperature of about 10'' K after 0.3 s as a result of hydrogen radiative 
cooling (see panel 2 of Fig. OJ. The initial beam impact causes the pressure to increase and a strong wave begins to 
propagate upward. The dynamics proceed slowly for the next three seconds, as much of the incident beam energy is 
radiated away. By 3.0 s, a significant fraction of the hydrogen in the region of beam energy deposition has been ionized 
and is no longer capable of radiating the beam energy, and so the temperature again quickly rises (see the hydrogen 
ionization fraction in panel 11 of Fig. 13). Over the next several seconds, the increasing temperature causes He i in 
the atmospheric range of 1.0 to 1.4 Mm to become ionized to He 11. The cooling increases through He recombination 
and the He 11 A304 line, and the radiative cooling again balances the beam and XEUV heating. When this balance is 
reached, a relatively gentle phase of the evolution begins. Meanwhile, the wave which resulted from the initial impact 
of the beam has been traveling through the atmosphere, bringing material from the chromosphere into the transition 
region and corona. By 50 s the initial wave has reached the top of the flux tube. Material pushed by the wave has 
increased the apex density by a factor of three and raised the temperature to about 2 MK. The increased coronal 
density causes the transition region to become much more extended in height than in the preflare atmosphere (panel 4 
of Figure In addition, with elevated coronal density the atmosphere becomes more effective at stopping the beam, 
which therefore deposits its energy progressively higher in the atmosphere. Column 4 of Figure 01 (panels 4, 8, 12) 
illustrates these effects at a representative time during the gentle phase. 

After FIO has evolved for about 73 s, most of the He 11 in the range of 2.0 to 2.3 Mm has been ionized, and the 
radiative cooling again decreases below the flare heating. This is illustrated in the first column of Figure 0] An 
explosive increase in temperature results and creates a supersonic wave which pushes material away from its front 
and forms a lower density (though still much higher than the prefiare density) high temperature bubble. The second 
column of Figure 0] (panels 2, 6, 10) shows these effects at 85 s. Meanwhile, more material has been evaporated into 
the corona and the loop apex density has increased approximately an order of magnitude over the preflare density 
(panel 10). For the next 150 s the bubble expands as the explosive wave propagates through the atmosphere, until the 
wave reaches the boundary of the fiux tube. The atmosphere attains a steady state, with a hotter, denser corona, and 
a transition region that is closer in height to the photosphere than in the preflare atmosphere (panels 4, 8, 12). The 
final apex temperature is 5.2 MK, the electron density is 5.9 x 10^° cm~^, and the mass density has increased by a 
factor of 70 reaching 1.2 x 10"^'^ g cm~"^. This temperature and density are somewhat lower than coronal values often 
observed for M-class flares. Although the XEUV heating rate would be larger from an even hotter, denser corona, it 
would still be much smaller than the beam heating (see below), and therefore our predictions of the optical and UV 
emission from the lower atmosphere are not affected. For example, an atmosphere with an apex temperature of 20 
MK and electron density of 5 x lO"'^ cm~'^ gives an XEUV heating rate which is five times larger than ours, but still 
50 times smaller than the beam heating rate. 

A major difference between the models presented here and those of AH99 is the time at which the explosive phase 
begins. The explosive phase of the FIO model of AH99 begins after about 27 s of impulsive heating compared to 73 
s for our FIO model. The reason is that the lower spectral indices we employ result in broader electron beam energy 
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3. FLARE SIMULATIONS 



3.1. FIO Flare Dynamics 
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deposition. With less concentrated heating, the atmosphere takes longer to ionize resulting in a longer delay before 
the onset of the explosive phase. 

3.2. Fll Flare Dynamics 

The initial evolution of the Fll flare is similar to FIO but proceeds much faster. Figure |S1 shows the atmosphere at 
several times during this flare. The beam rapidly heats the atmosphere in the region of impact to about 13,000 K, but 
plateaus after 0.04 s because hydrogen radiative cooling balances the beam heating. As in the FIO, case a high velocity 
wave begins to propagate upward from the initial beam impact. By 0.18 s, most of the hydrogen is ionized and the 
atmosphere again rapidly heats, but plateaus at just over 50,000 K as a result of He ii radiative cooling (first column 
of Fig.[SJ. By 1.0 s, helium is completely ionized and the radiative cooling is insufficient to balance the beam heating, 
causing the atmosphere to explosively heat (second column of Fig. EJ. A high-temperature, low-density bubble forms 
in a fashion similar to the FIO flare, but at a location much lower in the atmosphere than in the FIO case. Since the 
atmosphere heats so much more quickly in the Fll flare, the initial wave resulting from beam impact has not had 
time to evaporate much material into the corona. Therefore, the beam is still depositing most of its energy in the 
lower chromosphere when the explosive phase starts. The explosive bubble forms at about 1 Mm (panel 2 of Fig. EJ. 
A supersonic wave forms as a result of the explosion, and it is this wave which is primarily responsible for depositing 
chromospheric material into the corona, in contrast to the FIO case where the material was carried primarily by the 
initial wave. It is likely that for flares larger than Fll, the onset of the explosive phase will be extremely fast (i.e. 
hydrogen and helium will be ionized almost immediately). For very large flares, the explosive phase will essentially 
begin simultaneously with the beam heating, with no observable lag. 

Hawlev & Fisher (1994) speculated that XEUV emission could provide significant heating below the temperature 
minimum region, which might explain the increased continuum emission produced in white light flares. However, in 
both FIO and Fll, we find that the XEUV heating does not penetrate more deeply than the electron beam, nor does 
it contribute significantly to the energy deposition rate. Despite an order of magnitude more heating than in AH99, 
XEUV heating accounts for less than 1% of the total energy deposited in the FIO flare and about 7% in the Fll flare. 

3.3. Line Emission 

As material evaporates from the chromosphere into higher temperature regions, many emission lines brighten dra- 
matically. Light curves showing these increases in flux during the FIO flare are plotted in Figure |H| Lya quickly 
brightens as lower chromospheric material is heated to ~25,000 K. By 30 s, much of the plasma in the lower transition 
region has been further heated and Lya emission slowly begins to decrease. The He ii A304 line continually grows 
from the beginning of the flare as progressively more material is heated to its temperature of formation (^50,000 K) 
until the onset of the explosive phase. After the explosive phase begins, the expansion of the high temperature bubble 
ionizes most of the He II, and the intensity of He II A304 decreases (see panels 2 and 3 of Figure 0J. At the onset of 
the explosive phase a downward moving "condensation" wave is created and raises the plasma density in the lower 
transition region resulting in another increase in Lya emission. The peak occurring at 112 s is due to emission from the 
high density, lower temperature boundary of the expanding bubble (see panel 3 of Figure 0)), and the sharp decrease 
in intensity occurring at 185 s is the result of the bubble passing through the boundary of the flux tube. Ca ii K does 
not respond to the explosive phase until the condensation begins to increase the chromospheric density at about 115 
s. At 166 s, a large, downward-moving wave passes into the lower chromosphere increasing the density and causing 
Ca II K to brighten (compare the location of the transition region in panels 3 and 4 of Figure 01 . 

Figure shows light curves for the same emission lines during the Fll flare. As in FIO, the lines rapidly brighten in 
response to the flare heating. After 0.5 s, the region of beam impact has been heated past 40,000 K, and the rate of 
intensity increase slows for Lya and Ha. At 2.5 s, the explosive bubble begins to move to areas of lower density and 
the Lya and He ii A304 intensities slowly decrease. At 5.5 s, a condensation front passes into the lower chromosphere, 
increasing the density and causing Ca ii K to brighten. 

Mass motions during the flare cause significant Doppler shifts and introduce large asymmetries in the line profiles. 
Figure IHl shows profiles for the same emission lines, at successive times during FIO. In panels 2-4, the red wing of Lya 
is strongly enhanced. Panels 7-9 show Ha to be weakly red shifted (average Doppler shift of 2 km s^^). These indicate 
chromospheric condensation (mass motion downward toward the photosphere). However, the line center of He ii A304 
shows the opposite effect (central reversal in panel 13), with a blue shift exhibiting chromospheric evaporation. The 
line center of He ii A304 is formed higher in the atmosphere — where the velocity field is stronger — than the other 
lines shown and exhibits a peak upward v elocity of ^--^120 km s~^. A similar bl ue shift has be en observed in several 
solar flares l|Silva et alJll997tlBrosiu£ll2003l and references therein). In particular . iBrosiu^ l)2003j) observed a velocity of 
~140 km s~^ early in the impulsive phase of a moderate flare using an O ill line which forms at a similar temperature 
to He II A304. 

Figure |51 shows similar plots for line emission during the Fll flare. At 0.3 s both Ha and Ca ii K have switched from 
absorption to emission profiles, but neither Lya nor He ii A304 has significantly brightened. However, by 1.0 s the 
intensity of all these lines has markedly increased. In panel 10, He ii A304 has a blue shift corresponding to a velocity 
of ~ 60 km s-^ At 6.0 s there are red shifts present in Ca ii K and Ha indicating motion towards the photosphere. 
The region of Ha formation reache s a pea k downward velocity of about 40 km s~^ which is consistent with velocities 
reported bv llchimoto fc Kiirokawal ((1584!) ■ The profiles of Lya and He ii A304 are very asymmetrical since they are 
formed in regions with both strong upward and downward directed waves. It is, therefore, not possible to determine a 
single velocity from these profiles. In the last column, Lya and He ii A304 show enhanced blue wings, but the Ca ii K 
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line is red-shifted indicating a downward-moving wave in the lower chromosphere and an upward-directed wave in the 
upper chromosphere and transition region. 

Clearly, the detailed line profiles in the simulations d epend strongly on the fla re dynamics and evolve significantly 
during the flare. To interpret the line profiles, we follow ICarlsson fc SteinI l)1994fl and write the formal solution of the 
transfer equation for the emergent intensity as 

7,(0) = - / S^e-^-^l^dT, = - / S,^T,e-^"'^'dz (7) 

where Xi/ is the linear extinction coefficient, defined as the product of the density of emitting particles with their 
cross section. The integrand in this equation is known as the contribution function and indicates how much emission 
originates from a height, z. We divide the contribution function into three physically meaningful terms: Xvl'^v is large 
where there are many emitting particles at small optical depth, which is a situation that often arises in the presence of 
strong velocity gradients. 5*^ is the source function, defined as the ratio of the emissivity to the extinction coefficient. 
For a given line Sy is independent of frequency due to the assumption of complete redistribution. The last term is 
Tye^'^"!^ which represents the attenuation caused by optical depth and peaks at = 1. 

As an example, we consider the contribution function for the He ii A304 line at 73 s during the FIO flare. This is 
plotted in Figure [TUI in inverse gray scale (darker shades means higher intensity). The line profile, shown in the bottom 
right panel, has two large peaks in the wings on either side of line center. The upper right panel {Tye~'^''^^) shows 
that, in the wings, Ti, = 1 in the range of 1.1 — 4 Mm. In this range the source function and Xvl'^v are both large 
leading to bright emission in the wings from this atmospheric range. At line center, however, Ti, = 1 in a relatively 
small region around 5.2 Mm. At this height the source function is small resulting in a smaller contribution function 
and less emission. 

The He ii A304 line is formed in the transition region and is made complicated by the strong velocity gradients 
occurring there. In the lower chromosphere the velocity field is much smaller. As an example of emission resulting 
from this region, we consider, in Figure ITTl the contribution function for Ca ii K at 130 s after the start of the FIO 
flare. For line center, = 1 at a height of about 1.5 Mm. In this region Xvl^v and are also at their maximums 
indicating that most of the emission originates from this region. The wings are formed in the region of 0.3 — 0.5 Mm 
where Ti,e~'^"l ^ is at a maximum. However, the source function and Xvl'^v are smaller than at line center, so there is 
less emission in the wings. 

The line asymmetry of Ca ii K is more pronounced in the Fll flare (compare the last rows of Figure |H| and EJl . 
Figure [T^ shows the Ca ii K contribution function at 11.0 s after the start of Fll. Most of the emission originates from 
an area with a strong downward directed condensation wave. Comparison of this figure with Figure ITTl reveals that 
Ca II K is formed deeper and in a narrower region in Fll than in FIO, with a majority of the emission coming from a 
height of 0.78 Mm. During the Fll explosive phase the beam energy is deposited at lower height than in FIO. This 
has the effect of pushing the transition region to lower heights and compressing the chromosphere relative to FIO. 

3.4. Continuum Emission 

Observations of solar white light fla res show enhanced optical continuum emission that is temporally correlated 
with the impulsive phase of the flare ijHudson et al. 1 11992(1 . Our flare simulations also show dramatically increased 
contin uum emission includ ing the enhanced Balmer and Paschen jumps which are ubiquitous in type I white light 
flares fMetcal f et an i2003t). Figures 1131 and 1141 show continuum light curves for four spectral regions during FIO 
and Fll respectively. The Lyman continuum is very enhanced, reaching thousands of t imes its preflare valu e. This 
compares well with the increased Lyman continuum observed in an X-class solar flare bv iLemaire bmce 
the preflare intensity is much higher in the optical region, the relative increase of optical emission is moderate. At 5000 
A the optical continuum reaches a peak increase of about 5% for FIO and 30% for Fll. This result is very similar to 
observations. i Chen fc Dind l|2005j) observed, for a white light flare with a beam strength of 1.0 x 10^*^ ergs cm~^ s~^, a 
continuum enhancement of ~ 4 - 5%. In our simulations, both the Balmer and Paschen continua initially decrease in 
intensity (see the insets in Fig. ll3|l . This is discussed in AH99, and is due to the non-thermal ionization resulting from 
the electron beam. The beam raises collisional rates in the upper chromosphere causing a higher population density of 
excited states of hydrogen. This increases the probability that Balmer and higher order hydrogen continuum photons, 
which previously would have escaped, are absorbed, and results in a decrease in the continuum intensity. 

It is interesting to compare the total fluxes in the lines and continua, as shown in Figure [151 The beam and XEUV 
heating is also indicated. Despite differences in dynamics and duration, in both the FIO and Fll flares the continuum 
carries 78% of the total radiation and the lines carry the remaining 22%. 

3.5. Comparison to the AH99 Models 

As has been noted, these models differ from those of AH99 by including more realistic electron beam heating and 
greatly increased XEUV heating. This leads to a number of observational differences. Our electron beam penetrates 
deeper (see Fig. ^ and more directly heats the lower atmosphere. This results in a larger emitting region for the 
hydrogen Balmer lines (compare Fig. 13 with Fig. 3 of AH99) and more intense Balmer emission. This is especially 
important for the Balmer continuum, which is responsible for radiative backwarming below the temperature minimum 
region, where white light continuum emission is produced |Metcalf et al.. 1990) . Since our models have more intense 
Balmer continuum, they produce more backwarming, a hotter flaring photosphere, and more intense white light 
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emission than the models of AH99. For the Fll flare, AH99 found a 4% increase in white light emission, compared 
to the 30% increase we obtain in our Fll model. Another difference between our models and AH99 is the amount 
of He II A304 emission. Our spatially extended beam deposition produces a pronounced helium-emitting plateau in 
the transition region during Fll, and results in an order of magnitude brightening of the He ii A304 line compared to 
AH99. 

Despite the brighter emission in our simulations, the velocities and many line profiles are similar to those found in 
AH99. For example, in the Fll flare AH99 obtain a downward chromospheric condensation velocity of ~ 40 km s^^, 
which is comparable to our result. Also, the Ca II K line profile at 73 s in our FIO flare (panel 18 of Fig.lSJ is very 
similar to the profile calculated at 30 s in the FIO model of AH99 (panel 6 of Fig. 14 of AH99). As noted earlier our 
simulations evolve more slowly than those of AH99, so these line profiles occur near the beginning of the explosive 
phase in each simulation. 

4. LIMITATIONS AND PLANS FOR FUTURE MODELS 

We have created detailed models of the radiation and dynamics occurring during a moderate and a strong solar flare 
and have taken considerable effort to make them as accurate as possible. However, there are a number of limitations 
to these models. The atmosphere is assumed to be one-dimensional, plane-parallel and aligned with the magnetic 
field. In assuming a one-dimensional atmosphere, the possibility that radiation can escape the flux tube through the 
sides is neglected. The assumption of a circular loop geometry neglects the more tapered structure of true magnetic 
field loops and will affect the calculation of the thermal XEUV heating. Ions and electrons are assumed to have the 
same temperature. We also have assumed complete redistribution in calculating all transitions, and we have implicitly 
assumed that the gas pressure is lower than the magnetic pressure so that the gas is confined to the flux tube and 
MHD effects can be neglected. 

We are in the process of making a number of improvements to the models. The models shown in this paper represent 
a "general" flare, i.e. the electron beam parameters are assumed to remain constant throughout the flare. RHESSI 
hard X-ray spectral observations have sufficient time resolution to see temporal variations in the electron beam on the 
order of a few seconds allowing the evolving electron beam parameters to be inferred. Thus, simulations of specific 
flares may be constructed, and the calculated emission profiles can be directly compared to observations of lines and 
continua for that specific flare. A forthcoming paper will present results for a simulation of the large X-class 23 July 
2002 solar flare. 

5. CONCLUSIONS 

We have constructed radiative hydrodynamic simulations of the flaring solar atmosphere and have used them to 
explore the flare-induced optical and UV emission. These models extend the work presented in AH99 by using more 
realistic electron beam deposition rates and by incorporating backwarming from a large number of XEUV lines. As in 
AH99, we find that the impulsive fiare naturally divides into two phases, an initial gentle phase followed by a period of 
explosive increases in temperature and pressure. The explosive wave front creates a high temperature "bubble." The 
bubble expands until the wave passes through the boundary of the flux tube, at which point the atmosphere attains 
a steady state with a hotter corona and deeper transition region than in the preflare atmosphere. 

Both moderate and strong flares show large increases in line and continuum emission. As chromospheric material is 
heated to transition region temperatures, the He Ii A304 line is especially enhanced reaching values of many thousands 
of times its preflare level. The optical continuum also brightens showing a peak enhancement of 5% for FIO (comparable 
to observed enhancements) and 30% for Fll. The continuum dominates the emitted radiation, carrying 78% of the total 
emission for both FIO and Fll. Chromospheric evaporation causes blue-shifts in transition region lines corresponding 
to a peak velocity of ~ 120 km s^^. Chromospheric condensation increases the density in the chromosphere causing 
Ca II K and Ha to significantly brighten and resulting in an Ha red shift of ~ 40 km s~^. 

We incorporated the large number of transitions in the ATOMDB database to more accurately model XEUV back- 
warming of the lower atmosphere. We find that the XEUV heating is an order of magnitude larger than the soft 
X-ray heating of AH99. Nevertheless, the level of heating is still too small to significantly contribute to the flare 
dynamics. By comparing these models with those in AH99 we find that heating from a more realistic electron energy 
distribution leads to slower flare evolution, and more intense brightening of emission lines and continua. The larger 
Balmer continuum present in our simulations is especially important because it radiatively backwarms the lower solar 
atmosphere, enhanced white light emission. 

This work has been partially funded by NSF grant AST02-05875 and HST grants AR-10312 and GO-8613. The 
computations presented here were carried out on the Astronomy Condor Network at the University of Washington, 
and we would especially like to thank John Bochanski for contributing many hours of computer time. 
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TABLE 1 

Bound-Bound Transitions 



Atom Ajj (A) Transition Atom Ajj (A) Transition 



1215.70 


Lya 


Ca II 


3933.65 


2s - 


-2p3 


1025.75 






8662.16 


2d3 


-2pi 


972.56 


Ly7 




8498.01 


2dz 


-2P3 


949.77 


Ly5 




8542.05 


2d5 


- 2p3 


6562.96 


Ha 


He I 


625.58 


ls2 i5o - 


- Is 2s 


4861.50 


H/3 




601.42 


ls2 i5o - 


- Is 2s i5o 


4340.62 


H7 




10830.29 


Is 2s '^.Si 


- Is 2p '"'P^ 


18752.27 


Pa 




584.35 


Is^ KSo - 


Is 2p 


12818.86 


P/3 




20580.82 


Is 2s 


- Is 2p ipf 


40513.47 


Ba 


He II 


303.79 


Is 251/2 - 


Is 2p 251/2 


3968.46 


2s - 2pi 




303.78 


Is 25i/2 


- 2p 2p0 



TABLE 2 
Bound-Free Transitions 



Atom Ajc (A) Initial State Atom Aic (A) Initial State 



H I 911.12 n=l Ca II 1421.04 4p^-P3/2 

H I 3635.67 n=2 He I 503.98 ls2 ^So 

H I 8151.31 n=3 He i 2592.02 Is 2s ^Si 

H I 14419.07 n=4 He i 3109.80 Is 2s ^So 

H I 22386.68 n=5 He i 3407.63 Is 2p ^P9 

Ca II 1044.00 4s25^^2 He I 3663.37 Is 2p ipf 

Ca II 1217.50 3(^2 D3/2 He 11 227.84 Is 

Ca II 1218.40 3d2D5/2 He 11 911.34 2s 251/2 

Ca II 1416.55 4p2p0^^ He 11 911.36 2p 2p0^^ 



TABLE 3 
Soft X-ray EUV Wavelength Bins 



Wavelength Range (A ) 


Central Wavelength (A ) 


1 - 15 


11.2 


15 - 25 


18.6 


25 - 50 


37.2 


50 - 100 


72.3 


100 - 150 


125.5 


150 - 200 


177.8 


200 - 250 


224.7 


250 - 300 


269.2 


300 - 500 


352.4 


500 - 750 


608.5 


750 - 1000 


867.7 


1000 - 1500 


1163.6 


1500 - 2000 


1695.4 


2000 - 2500 


2362.1 
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Fig. 1. — Electron beam heating rates in the prefiare atmosphere for the FIO flare. The soUd hne shows the heating rate used in this 
paper (defined in Eq. |^ and the dashed line sh ows t he initial beam heating rate of AH99. Our use of recently measured spectral indices 
and cutofi^ energy for the beam as described in S 12.21 results in extended beam heating that penetrates deeper in the atmosphere than the 
models described in AH99. 
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Fig. 2. — The XEU V heating rate as a function of atmospheric height compared with the soft X-ray heating rate of AH99 at a comparable 
time in the late stages of a moderate flare. 
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FIO at t= 0.04 s F10att=0.3s F10att=3.0s FIO at t= 50.0 s 




Height (Mm) Height (Mm) Height (Mm) Height (Mm) 

Fig. 3. — The solar atmosphere at four times during the gentle phase of the FIO flare. The top row shows the log of the temperature, 
T as a function of height compared with the preflarc state (PFS). In the middle row the electron density, ne (left axis) and beam heating 
rate, Qe (right axis) are plotted. The bottom row shows the mass density, p (left axis) and hydrogen ionization fraction, X, (right axis). 
Note the change in scale of the horizontal axis in the last column. 
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FIO at t= 73.0 s FIO at t= 85.0 s FIO at t= 150.0 s FIO at t= 221.0 s 




2468 10 2468 10 2468 10 2468 10 
Height (Mm) Height (Mm) Height (Mm) Height (Mm) 



Fig. 4. — The FIO solar atmosphere at four times during the explosive phase. The quantities plotted are identical to Fig. ISlexcept that 
the He III fraction is plotted in the bottom row rather the hydrogen ionization fraction. 



14 



AUred et al. 



Fllatt=0.6s Fllatt=1.0s Fllatt=6.0s Fll at t= 15.2 s 




Height (Mm) Height (Mm) Height (Mm) Height (Mm) 

Fig. 5. — The solar atmosphere at four times during the Fll flare. The quantities plotted are identical to Fig. 1^ except that the He III 
fraction is plotted in the bottom row rather the hydrogen ionization fraction. 




50 100 150 200 

Time (s) 



Fig. 6. — Light curves for Lyo, Ha, He II A 304 and Ca II K during the FIO flare. The preflare emission has been subtracted. The line 
flux for Ca II K includes only the central reversal. The explosive phase begins at 73 s and is marked with a dashed line. 
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2 4 6 8 10 12 14 

Time (s) 



Fig. 7. — Light curves for Lyo, Ha, He II A 304 and Ca II K during the Fll flare. The notation is identical to Fig. |S] The explosive 
phase begins at 1 s and is marked with a dashed line. 
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Wavelength from Line Center (A) 

Fig. 8. — Line profiles for Lya, Ha, He II A 304 and Ca II K at five times during the FIO flare (times are indicated at the top of each 
column). The dotted lines indicate the level of the continuum and line center, while the dashed line is the preflare line profile. The emission 
in panels 1 and 11 is too small to be seen on this scale. 
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o 

Wavelength from Line Center (A) 

Fig. 9. — Line profiles for Lya, Ha, He II A 304 and Ca II K at four times during the Fll flare (times are indicated at the top of each 
column). The dotted lines indicate tlic level of the continuum and line center, while the dashed line is the preflare line profile. The emission 
in panels 1 and 9 is too small to be seen on this scale. 
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500 -500 500 -500 

Al/ (km/s) l\v (km/s) 



Fig. 10. — Components of the contribution function for the He li A304 line at 73.0 s from the start of the FIG flare, are plotted in inverse 
gray scale (high intensity is darker). The vertical axis is atmospheric height, and the horizontal axis is frequency from line center in units 
of km s~^. The components plotted are Xvl^v (upper left), T^e"''"/'' (upper right), 5„ (lower left), and the contribution function, which 
is the product of the previous three (lower right). For reference, the line profile has been scaled and plotted in the lower-right panel. The 
source function and Planck function at line center are plotted in the lower left panel to show the extent to which the atmosphere diverges 
from LTE. In addition, lines indicating Ti/ = 1 are plotted in the top two and bottom left panels, and the plasma velocity is plotted in the 
top two and lower right panels. Positive velocity indicates material moving outward, i.e. in the direction of the corona and correspond to 
blue shifts in the line profile (see A scale on top axis). 
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Fig. 12. — Components to the contributi on f unction for the Ca II K Hne at 11.0 s from the start of the Fll flare. The plot is in inverse 
gray scale. The notation is identical to Fig. 1101 
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50 100 150 200 250 

Time (s) 

Fig. 13. — Light curves are shown for four continuum regions during the FID flare. The preflare continuum has been subtracted. The 
insets in the bottom two panels show the first second of the flare, and illustrate the initial continuum dimming. The major contributors 
to emission in each bin arc He I and He II continua in the top bin; Lyman and He II continua in the second bin; Balmer, He I, and Ca II 
continua in the third bin; and Paschen and He I continua in the bottom bin. 
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Fig. 14. — Light curves are shown for four continuum regions during the Fll flare. The notation is identical to Fig. 1131 Insets show the 
first 0.5 s of the flare and illustrate the initial continuum dimming. 
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Fig. 15. — Light curves of the total line and continuum fluxes are given for FIO (a) and Fll (b). The dotted and dashed lines indicate 
totals over continua and lines respectively, and the solid line is the sum of lines and continua. In each case, the preflare emission has been 
subtracted so that the plots represent the emission resulting from the flare heating. The dot-dashed line indicates the level of non-thermal 
beam heating, and the dashed-dot-dot line in panel (b) is the XEUV heating for Fll. For the FIO flare, the XEUV heating is too small to 
be seen on this scale. 



